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ABSTRACT 

We present multi-wavelength observations of the centre of RXCJ1504. 1-0248 - the 
galaxy cluster with the most luminous and relatively nearby cool core at z ^ 0.2. 
Although there are several galaxies within 100 kpc of the cluster core, only the brightest 
cluster galaxy (BCG), which lies at the peak of the X-ray emission, has blue colours 
and strong line-emission. Approximately 80MQyr~^ of intracluster gas is cooling 
below X-ray emitting temperatures, similar to the observed UV star formation rate of 
~140 M0 yr~^ . Most star formation occurs in the core of the BCG and in a 42 kpc long 
filament of blue continuum, line emission, and X-ray emission, that extends southwest 
of the galaxy. 

The surrounding filamentary nebula is the most luminous around any observed 
BCG. The number of ionizing stars in the BCG is barely sufficient to ionize and heat 
the nebula, and the line ratios indicate an additional heat source is needed. This heat 
source can contribute to the Ha-dcduccd star formation rates (SFRs) in BCGs and 
therefore the derived SFRs should only be considered upper limits. 

AGN feedback can slow down the cooling flow to the observed mass deposition 
rate if the black hole accretion rate is of the order of 0.5 M0 yr"-'^ at 10% energy 
output efficiency. 

The average turbulent velocity of the nebula is Wturb ~ 325 kms^^ which, if shared 
by the hot gas, limits the ratio of turbulent to thermal energy of the intracluster 
medium to less than 6 percent. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: RXCJ1504.1- 
0248 - cooling flows - ionized gas 



1 INTRODUCTION 

Cool core clusters are characterised by bright central peaks 
in their X-ray surface brightness profiles, cool and radially- 
increasing core temperatures, and short central cooling 
times. Unless the gas is reheated, the large amount of en- 
ergy radiated in X-rays implies several hundred solar masses 
of gas cools below lO'' K every year (Cowie & Binney 1977; 
Fabian & Nulsen 1977). However, X-ray spectra from XMM- 
Newton and Chandra show cool-core clusters lack strong 
cooling lines, so less gas is cooling below X-ray emitting 
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temperatures than is expected from the X-ray luminosity 
(Peterson et al. 2001, 2003; Bohringer et al. 2001; Peterson 
& Fabian 2006; McNamara & Nulsen 2007). The implication 
is that the intracluster medium (ICM) must be heated. 

On the other hand, large quantities of cool gas and dust 
have been found in brightest cluster galaxies at the centre 
of cool-core clusters. Searches conducted at optical and in- 
frared wavelengths reveal both warm and cool gas reservoirs 
(e.g. Edge 2001; Salome & Combes 2003; Hatch et al. 2005; 
Jaffe, Bremer & Baker 2005; Johnstone et al. 2007). The 
BCGs in cool-core clusters generally have anomalous blue 
colours implying recent star formation, and emit in recom- 
bination and low-ionization lines (McNamara 1997; Peres 
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et al. 1998; Cavagnolo et al. 2008). These observations im- 
ply a relationship between the properties of the cluster and 
the central galaxy, but this interaction has not yet been un- 
derstood or quantified precisely. 

RXCJ1504. 1-0248 is one of the most massive cool 
core clusters, with a prominent central X-ray brightness 
peak, and a short central cooling time. A classical cool- 
ing flow model leads to a mass deposition rate of 1400 — 
1900 M© yr~^ , and the brightest cluster galaxy emits strong 
low ionization emission lines (Bohringer et al. 2005). These 
properties mark this cluster as a prime target to study the 
relationship between the ICM, the brightest cluster galaxy, 
and a possible central active galactic nucleus (AGN). 

Here we report results of a multi-wavelength study of 
the core of RXCJ1504. 1-0248. It lies at z = 0.2153, which al- 
lows for a detailed spatial analysis of the BCG. We present 
results of integral field unit (IFU) spectroscopic observa- 
tions aimed at understanding the morphology, kinematics 
and ionization state of the line-emitting nebula surrounding 
the BCG. We also present an analysis of X-ray spectra ob- 
tained with the XMM-Newton Reflection Grating Spectrom- 
eter (RGS) and European Photon Imaging Cameras (EPIC) 
and compare the ICM mass deposition rate to the star for- 
mation rate of the BCG. Sections 2 and 3 describe the ob- 
servations and the data reduction process, while in section 
4 we present the properties of the BCG, the nebula, and 
the cluster core cooling rate. In section 5 we combine the 
results from the multi-wavelength observations to quantify 
the relationship between the ionized gas, the BCG and the 
cluster core. We summarize our findings in section 6. 

Unless otherwise stated, the following cosmological pa- 
rameters have been used: Q.m = 0.3, SIa = 0.7, Ho = 
70kms~^Mpc"\ For a redshift z = 0.2173 (the redshift 
of the BCG), the scale is 3.519 kpcarcsec"^. 



2 OBSERVATIONS 

2.1 Optical images 

We observed RXCJ1504. 1-0248 on 2006 March 1 with the 
Visible Multi-Object Spectrograph (VIMOS; LeFevre et al. 
2003) on the 8.2-m UTS of the Very Large Telescope (VLT) 
at ESQ Paranal in Chile, operated in imaging mode, using 
the R-band filter. The VIMOS field of view consists of four 
optical channels (also referred to as "quadrants"), which in 
imaging mode cover a field of view of approximately 7' x 8' . 
The CCD size is 2048 x 2440 px^, therefore yielding a pixel 
size of 0.205 arcsec. Twelve scientific exposures were taken, 
each with an integration time of ~ 90 s. 

2.2 Optical integral field spectroscopy 

On 2006 April 5-6, the central galaxy of the cluster was also 
observed with VIMOS in Integral Field Unit (IFU) mode 
with the LR-red and LR-blue grisms. Both grisms provide 
a 54 X 54 — arcsec^ field of view, divided into 6400 optical 
fibres that offer a spatial sampling of 0.67 arcsec/fibre. The 
fibre to fibre distance on the sky is about 1 arcsec and the 
dead space between fibres is below ~0.1 arcsec. Thus, the 
optical fibres, coupled to an array of 6400 microlenses, en- 
sure a nearly continuous sky coverage. All quadrants of the 



Table 1. Summary of the VIMOS IFU observations. 

Wavelength Exposure Seeing 

Grism Filter coverage (A) time (sec) (aresee) 

LR-red OS-red 5900 - 9150 8 x 450 0.5 - 1.3 

LR-blue OS-blue 4000 - 6700 8 x 450 0.5 - 0.8 

spectrograph are divided into four pseudoslits, each orga- 
nized into a 20 x 20 array of microlenses and holding five 
20 X 4 fibre modules. 

The LR-red grism, used together with the OS-red filter, 
samples the wavelength range 5900 — 9150 A, with a dis- 
persion of approximately 7.3 A pixel" ^ and a spectral res- 
olution of 260. The LR-blue grism, used together with the 
OS-blue fllter, covers the wavelength range 4000 — 6700 A, 
with a dispersion and spectral resolution of approximately 
5.3 ApixeP^ and 220, respectively. Eight 450-second expo- 
sures were taken in each setting, at averaged airmasses in 
the range 1.093 — 1.717, and with a pointing dither of about 
(Aa, A(5) = (4-3 arcsec, — 3 arcsec). A summary of these ob- 
servations is provided in Table 1. 

2.3 X-ray observations 

RXCJ1504. 1-0248 was furthermore observed with XMM- 
Newton on 2007 January 22, for 39 ks. We extracted a 
lightcurve for each of the EPIC detectors separately and ex- 
cluded the time periods in the observation when the count 
rate deviated from the mean by more than 3(t in order to 
remove flaring from soft protons (Pratt & Arnaud 2002). We 
did not find significant time intervals affected by flares, such 
that after this cleaning, the net effective exposure was ~37.4 
ks and ~ 28.5 ks for EPIC/MOS and EPIC/pn, respectively. 



3 DATA REDUCTION 

3.1 VIMOS Imaging 

The imaging frames were reduced with the VIMOS Pipeline 
Recipes package, version 1.0, following all the steps de- 
scribed in section 5.5 of the VIMOS Pipeline User's Guide^, 
including flat fielding and bias subtraction. Flux calibration 
was applied using observations of standard stars in the L98 
field^ taken on the same night as the science frames. 

3.2 VIMOS IFU 

3.2.1 Reduction of the IFU data cubes 

The data reduction process was performed using the VI- 
MOS Interactive Pipeline and Graphical Interface version 

1.3 (Scodeggio 2005). Both the red and the blue datacubes 
underwent the usual VIMOS IFU reduction steps (Zanichelli 
2005; Covone et al. 2006), including adjustment of the "first 
guess" parameters of the instrumental model in order to cor- 
rect for optical and spectral distortions, spectra extraction, 
cosmic-ray cleaning, wavelength calibration and flux calibra- 
tion. The LR-red data were calibrated using the standard 

^ ftp:/ /ftp. eso.org/pub/cpl/vimos/vimos-pipeline-manual.pdf 
^ http:/ /wwwS.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/STETSON/standards/ 
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star EG-274, observed immediately after the science obser- 
vations were taken. Unfortunately, no standard star was ob- 
served at the same time as the observations taken with the 
LR-Blue grism. Therefore, the flux calibration of the LR- 
blue spectra was done using observations of the standard 
star LTT-3864 taken on 2006 May 27. 

Each fibre of the IFU has a slightly different transmis- 
sion. To obtain accurate fluxes we corrected for different 
transmission between fibers. We measured the flux in the 
7246 A sky line in the red data and 5577 A sky line in the 
blue data, and the spectrum from each flbre was normalized 
so that the flux in the sky line equals the median of the sky 
line flux from all the fibres. At this stage, spectra from fibres 
with low transmission were removed as they were too noisy. 

A sky spectrum was created for each fibre from the me- 
dian of 20 nearby spectra in the red frame (or 30 fibres in 
the blue frame) which contained little or no line emission. 
These sky spectra were then subtracted from the calibrated 
spectra. Finally the spectra were trimmed to the wavelength 
range 3885-6225 A and 5635-8348 A, for the blue and red 
gratings respectively. 

The absolute fiux calibration of the IFU data after the 
VIPGI reduction is uncertain as our total flux measurement 
of [On] is a factor of ^ 2 lower than the slit spectrum of 
(Bohringer et al. 2005). Therefore we fiux calibrate the IFU 
data using all stars and galaxies of known magnitude within 
the field-of-view We find the SDSS magnitudes of these ob- 
jects are consistent with those measured from our VIMOS 
image. The IFU spectra of the stars and galaxies are con- 
volved with g', r' and i' and VIMOS R-band filters, and the 
flux calibration is determined using SDSS magnitudes. This 
self-consistent calibration results in spectra which are 2.5 
times brighter than spectra calibrated using VIPGI. 

To increase the signal-to-noise ratio the spectra were 
binned using a Voronoi-2D binning algorithm'^ (Cappellari 
& Copin 2003). We use the weighted Voronoi tessellation as 
laid out by Diehl & Statler (2006) to have a signal-to-noise 
ratio of 15 in the H/3 line flux measured in the blue frame. 



3.2.2 Spectral emission line fitting 

The emission lines were measured by fitting a linear contin- 
uum plus Gaussian model profile to each line using the idl 
package MPFIT which performs a Levenberg-Marquardt 
least-squares fit. 

The doublet [On] AA3727, 3729 is never resolved, there- 
fore it is fit with a single Gaussian. The Ha line and the 
nearby [Nil] AA6548, 6583 doublet are resolved and the com- 
plex is fitted with three Gaussians, with the integrated flux 
of the [Nil] A6548 line fixed to be a third of the [Nil] A6584 
line. Similarly [Ol]A6363 is fixed to be a third of the fiux of 
[Ol]A6300, and [Olll]A5049 is fixed at a third of the flux of 
[Oiii]A5007. 

In the red wavelength domain, the spectra are strongly 
affected by fringing beyond 7000A. Consequently the signal- 
to-noise ratio at the wavelength of Ha is not as high as that 
achieved at the wavelength of [Oil]. 



^ using IDL codes from http: / /www- 

astro. physics. ox. ac.uk/'^mxc//idl/#binning 

* http://cow.physics.wisc.edu/~craigm/idl/fitting.html 
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Figure 1. Extinction map [E{B — V)] of RXJ 1504 derived from 
the ratio of Ho/H/3 measured in each bin. The centre of the galaxy 
is marked by a cross. Black bins indicate areas where [On] is 
detected, but H/3 is either not detected or the Ha/H/3 ratio implies 
E{B - V)=0. 



3.2.3 Extinction along the line of sight 

We assume that there are two sources of dust which lead to 
extinction and reddening of the spectra. The first source is a 
screen of foreground dust at the redshift of the EGG and the 
other source is dust within our Galaxy. We first correct the 
spectra for Galactic extinction of E(B-V)=0.108 (Schlegel 
et al. 1998) , then blue-shift the spectra to the rest-frame of 
the BGG, and correct for extinction caused by dust in the 
EGG. 

The extinction caused by dust within the BGG is deter- 
mined from the observed Ha/H/3 ratio. The Ha/H/3 fiux is a 
good indicator of extinction along the line of sight, as shorter 
wavelengths are scattered stronger by dust particles. There- 
fore, in dusty regions, values of Ha/H/3 are significantly 
above the theoretical ratio of 2.87 derived by Osterbrock & 
Ferland (2006) for Case B recombination at Tc = 10, 000 K. 
We assume the H/3 flux is depleted due to 2A of H/3 stel- 
lar absorption (the mean found in a sample of star forming 
galaxies, Buat et al. 2002) and we correct the Ha/H/3 ratio 
for this before calculating the extinction. 

The E{B — V) colour excess in each bin was com- 
puted using the equations from Moustakas, Kennicutt Jr. 
& Tremonti (2006), 
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where fcna — kua ~ 1.278 for Rv=3.2. The spatial distribu- 
tion of the extinction is plotted in Fig. 1. The central region 
of the galaxy is dustier than the outskirts, but the dust is 
patchy and there is no overall smooth pattern. The spec- 
tra in each bin was blue-shifted to the rest-frame and then 
corrected for extinction using the equation of Galzetti et al. 
(2000) with the Rv = 3.2 extinction law. 
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3.3 Reduction of the XMM-Newton data 

We reduced the XMM-Newton data using tlie 7.0.0 version 
of the Science Analysis System (SAS); the standard analysis 
methods using this software are described in e.g. Watson 
et al. (2001). 

We subtracted the instrumental background from the 
EPIC spectra using closed-filter observations which were 
normalized to match our observation using the count rates 
in the hard energy band (10 - 12 keV for MOS, 12 - 14 
keV for pn) outside of the field of view (OoFoV). Out-of- 
time events were subtracted from the EPIC/pn data using 
the standard SAS prescription for the extended full frame 
mode. The cosmic X-ray background (CXB) obtained from 
spectra extracted from an annulus between 9.5' and 12.5' 
(outside the cluster's rsoo) is modeled with three compo- 
nents: two thermal components to account for the local hot 
bubble (LHB) emission (fcTi = 0.08 keV) and for the Galac- 
tic halo (GH) emission {kT2 ~ 0.2 keV), as described by 
Kuntz & Snowden (2000), and a power-law to account for 
the integrated emission of unresolved point sources and for 
possible contamination from the residual soft proton par- 
ticle background. Since different detectors can be affected 
differently by soft protons, we leave the power-law indices 
and normalizations free between EPIC/MOSl, M0S2 and 
pn. The temperature of the LHB was frozen to 0.08 keV, 
while the temperature of the GH and the spectrum normal- 
izations of the thermal components were free in the fit, but 
constrained to be the same for all three EPIC detectors. The 
spectra were fit in the 0.4-7. keV energy band. 

The RGS spectra were extracted following the method 
described by Tamura et al. (2001). We modeled the back- 
ground using the standard background model available in 
SAS (rgsbkgmodel, Gonzalez-Riestra 2004). The cluster 
spectra were extracted from a region which is 2' wide in the 
cross-dispersion direction of the instrument. The line emis- 
sion observed with the RGS from extended sources is broad- 
ened by the spatial extent of the source along the dispersion 
direction. In order to account for the line broadening in the 
spectral modeling, we convolve the line spread function (Isf ) 
model with the surface brightness profile of the source along 
the dispersion direction derived in the 0.8-1.4 keV band. We 
fit the 1st order spectra in the wavelength band of 8-25 A 
and the second order spectra in a narrower band of 8-16 A. 
The 1st and 2nd order spectra obtained by the two RGS de- 
tectors were fit simultaneously with the relative instrument 
normalizations left as free parameters. 

Freezing the Galactic absorbing column density (nn) 
to 6 X 10^" cm~^, as determined from Dickey & Lockman 
(1990) from 21 cm observations, gives very poor fits to the 
cluster EPIC spectra. Therefore, we froze the nn to 8.7 x 
lO'^" cm~^, the best-fit value determined by fitting EPIC 
spectra of the central 5' in combination with ROSAT All- 
Sky Survey (RASS) spectra to better constrain the emission 
in the soft X-ray band. This value agrees within the 90% 
confidence interval with the nn value computed from the 
100 ^m IR data, using the Boulanger et al. (1996) IR-tih 
correlation function (more details will be provided by Moi§ 
et al., in prep). 



Table 2. XMM-Newton spectral fit results for tlie central region 
of RXCJ 1504.1-0248. The abundance values are given with re- 
spect to solar abundances (Lodders 2003). Errors are quoted at 
the la level. 
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kT (keV) 
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c R+0.06 


O/Fe 


0.47 ±0.25 


0.47 (fixed) 


Ne/Fe 
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0.71 (fixed) 


Mg/Fe 
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0.57 (fixed) 


Si/Fe 




0.79 ±0.12 


S/Fe 




0.53 ±0.15 


Fe 


-I 77+0.94 
^- "-0.56 


47+0-01 


M (M0/yr) 




78l?o 



/ d.o.f. 467/481 3003/2494 



4 RESULTS 

4.1 The mass deposition rate from X-ray spectra 

To determine the amount of ICM gas cooling out of the 
X-ray emitting temperature range, we extracted the EPIC 
spectrum from a circle of radius 140 kpc (0.67') centred on 
the cluster centre. This is the cooling radius obtained by 
Bohringer et al. (2005) for a Hubble constant Ho of 70 km 
s-i Mpc-\ 

The spectrum is fitted with a Galactically absorbed 
single temperature coUisionally ionized optically thin mekal 
plasma model plus a classical cooling fiow model (with the 
lower cutoff temperature fixed to 0.1 keV). To ensure a cor- 
rect propagation of errors due to uncertainties in background 
determination, the EPIC source spectra are fit in parallel 
with the CXB spectra obtained from an outer annulus. The 
normalizations of the cluster spectral components were set 
to zero for the CXB data sets, while the normalizations of 
the background models were fixed for the source analysis to 
a ratio corresponding to the relative sizes of the extraction 
regions for the source and CXB only. 

The metal abundances were coupled between the sin- 
gle temperature and cooling flow models both for fitting the 
EPIC and the RGS data. The O/Fe, Ne/Fe and Mg/Fe ra- 
tios were fixed in the EPIC fit based on the best-fit values 
obtained from the RGS spectra, since the energy resolution 
and effective area calibration of EPIC around the energy of 
the O, Mg and Ne lines are too poor to allow a reliable de- 
termination of their abundances with EPIC alone, especially 
for such a hot cluster. The Si and S abundances of the hot 
plasma can only be measured with EPIC because the emis- 
sion lines of these elements lie outside the RGS energy band. 
The fit results are summarized in Table 2. Errors are quoted 
at the la level. Abundance ratios are given in proto-solar 
units (Lodders 2003). 

The high luminosity and the highly peaked surface 
brightness distribution of the cooling core in RXCJ1504.1- 
0248 allow us to obtain relatively sensitive XMM-Newton 
RGS spectra of this cluster despite its large distance. We 
show in Fig. 2 the spectrum obtained by combining data 
from the two RGS detectors. This is one of the most distant 
galaxy cluster spectra obtained with RGS. The cause for the 
discrepancy between the temperatures measured with EPIC 
and RGS is the different spectral extraction region (the RGS 
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Figure 2. Combined RGS 1st and 2nd order spectrum of 
RXCJ1504. 1-0248 and best-fit model. 



extraction region is effectively lO'-long in the dispersion di- 
rection and 2'-wide in the cross-dispersion direction, and the 
temperature profile of the cluster rises steeply outside the 
cooling radius of 0.67'). Because of the higher temperature 
measured by RGS, the measured absolute abundances are 
also high, but the abundance ratios are typically robust. 

The best-fit values of the mass deposition rates from the 
EPIC and RGS data are 78 and 68 Mq yr"^ respectively. 
The 3cr upper limits are 156 and 239 Mq yr~^, respectively. 

The so-called "cooling flow problem" remains in that 
the mass deposition rates obtained from X-ray spectra are 
significantly less than those inferred from the X-ray luminos- 
ity of between 1400 and 1900 Mq yr"^ for Ho = 70 km s~^ 
Mpc~^. This implies that significant heat input is needed to 
reduce the cooling from 1400-1900 down to at most 150-240 
solar masses per year (~10%). The bolometric luminosity 
calculated by Bohringer et al. (2005) is 4.3 x 10*^ ergs s~^, 
out of which more than 70% {Lxc = 3 x 10*^ ergs s~^) is 
emitted from inside the cooling radius. Assuming 10% of the 
gas is cooling and forming stars, more than 90% of Lxc must 
thus be supplied by a heat source, amounting to 2.7 x 10*^ 



ergs s 



4.2 Galaxies within the cluster core 

The cluster core contains the BCG and a number of 
fainter galaxies marked B1-B6 in the VIMOS R-band im- 
age (Fig. 3). Spectra of these galaxies are shown in Fig. 3 
and the BCG spectrum in Fig. 4. The light from Bl fell on 
some bad IFU fibres so we were not able to extract a spec- 
trum. All spectra were corrected for Galactic extinction of 
E{B-V)= 0.108 (Schlegel et al. 1998) and the BCG was cor- 
rected for additional intrinsic extinction of E{B — V)= 0.211 
derived from the measured Ha/H/3 ratio of 3.68 (assum- 
ing 2 A of H/3 absorption). The BCG emits strongly in re- 
combination and low ionization lines, and the observed and 
extinction-corrected emission-line ratios are given in Table 
4. 

Redshifts of galaxies B2-B6 were measured using the 
strong stellar absorption lines Ca ll H and K near 4000A. 
These lines are not visible in the BCG spectrum so the [On] 
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Figure 4. Extinction-corrected composite spectrum of the ion- 
ized nebula surrounding the BCG. The full height of the [On] 
emission line is not shown. Emission line ratios are given in Ta- 
ble 4. 



emission line is used instead. Redshifts are listed in Table 3, 
together with the strength of the 4000A break (D4000) and 
the rest-frame B — R colour. 

Galaxies B2, B3, B5 and B6 lie within a few hundred 
kms^^ of the redshift of the BCG therefore these galaxies 
are situated within the cluster. B4 is blue-shifted by almost 
2000 kms^^ compared to the BCG and is either a galaxy 
in the infall region of the cluster or a foreground galaxy. 

D4000 is defined as the ratio of the average contin- 
uum between the rest-frame 4050-4250A and 3750-3950A 
(Bruzual A. 1983). However, this wavelength range includes 
strong emission lines visible in the spectrum of the BCG. To 
remove emission-line contamination from the D4000 mea- 
surement, we alter our definition for the BCG to the ratio 
between rest-frame 41 15-4250 A and 3750-3840A. 

The strength of this break, D4000, is typically greater 
than 1.6 in elliptical galaxies (Kauffmann 2003), even if 
the elliptical galaxies reside at redshifts of up to z 1.1 
(Pasquali et al. 2006). The cluster galaxies B2-B6 have 
D4000 values typical of ellipticals, but the BCG has a much 
lower D4000, similar to late-type galaxies. 

B2 is the only galaxy, apart from the BCG, to exhibit 
[On] line emission. Although it lies close to the emission- 
line nebula that surrounds the BCG, the kinematics of the 
emission line gas (see section 4.4.4) imply the line is emitted 
from gas within B2 and not from the BCG nebula. 



4.3 The brightest cluster galaxy 

4.3.1 Mass 

The BCG of RXCJ1504.1-0248 is a massive galaxy: its A'^ag 
is 13.130 (2MASS). The mass-to-light ratios of Cappellari 
et al. (2006) imply a stellar mass of ~ 6 x 10^^ Mq. This 
is twice the mass of the nearby cD galaxy M87. Although 
no K-correction has been applied, the galaxy is red at these 
wavelengths with H ~ K — 0.546, so the K-corrected -K^mag 
may be even larger. The dominant source of error in the 
derived mass is likely to be the large scatter in the mass-to- 
light ratio. 
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Figure 3. Left: i?,— band VIMOS image of tlie entire field-of-view of tlie VIMOS IFU centered on the brightest cluster galaxy (labelled 
BCG). A number of other nearby bright galaxies have been labelled B1-B6. Right: Spectra of galaxies B2-B6. Most absorption features 
above 6OOOA are due to atmospheric absorption. Stellar absorption Ca K and H lines are detected in all galaxies and [Oil] is seen in the 
spectrum of B2. 



Red continuum 



Blue continuum 



-18.5 -18.0 -17.5 -17.0 -16.4 -15.9 -15.4 




(b) 





10 











.s- 


-10 




-20 

-30 




[Oil] flux 



10 20 30 
kpc 

B-R colour 



10 20 30 40 
kpc 

D4000 



-17.0 -16.3 -15.7 -15.0 -14.3 -13.7 -13.0 -0.7 -0.6 -0.4 -0.2 -0.1 0.0 0.2 0.6 0.8 1.1 1.3 1.5 1.S 2.0 





-20 -10 10 20 30 40 
kpc 



-20 -10 10 20 30 40 
kpc 



-20 -10 10 20 30 40 
kpc 



Figure 5. Continuum images of the BCG from left to right: (a) unsharp-mask of the i?— band VIMOS image showing the NE and SW fil- 
ament; (b) continuum from the rest-frame ~5400A in log units of ergcm"^ A ^ arcsec^^; (c) continuum from the rest-frame ~3200A 
in log units of ergcm"^ s~^ A ^ arcsec"^; (d) map of the [On] emission with colourbar in log scale in units of ergcm"^ s~^ arcscc"^; 
(e) rest-frame B — R continuum colour and (f) D4000. The noisy regions (where the flux per binned region is less than the rms noise) 
are coloured black in (b)-(d) and white in regions (e) and (f). The galaxy B2 is located at (-20,0) in the IFU images. 
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Galaxy 


BCG 


B2 


B3 


B4 


B5 


B6 


Redshift 


0.2173 


0.2176 


0.2151 


0.210 


0.221 


0.216 


Velocity shift relative to BCG (kms^^ ) 




65 


-550 


-1885 


800 


-385 


Projected distance (kpc) 




24.6 


26.3 


70.7 


68.8 


85.5 


D4000 


1.28 (1.11)^ 


1.77 


3.84 


1.79 


2.85 


1.80 


B — R colour^ (mag) 


-0.12 (-0.40)2 


0.0 


0.13 


0.43 


0.43 


0.23 



Table 3. Properties of the brightest cluster galaxy and galaxies nearby in projection. Redshifts are determined from the Ca II H and K 
absorption lines. Projected distances are given in kpc, as observed from the dust-uncorrectcd continuum. B — R colours are measured from 
regions free from any emission line contamination. ^ The B — R colour is the mag nitudc difference between 5600-5800A and 6900-7100A. 
2 Values in the brackets are derived from the extinction-corrected spectrum. 
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Figure 6. Emission line ratios of the nebula of the BCG from left to right: [Olll]/H/3, [Oll]/H/3, [Nll]/Ha. Binned regions are only 
displayed where [Olll], H/3, [Nil] and Ha fluxes arc greater than 5x10^"*^^ ergcm^^s^-*^ arcsec^^. All scales are in log. 
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Figure 7. Optical spectrum of RXCJ1504. 1-0248 with emission- 
lines masked out (black solid line) and the best-fit models over- 
plotted in the dashed line. The top panel shows a 2-component fit 
to the SED consisting of an elliptical galaxy template plus a B3V 
stellar template. The coloured lines show the relative luminosity 
from each component. The bottom panel shows a 3-component fit 
to the SED consisiting of an elliptical galaxy template, an 05V 
stellar template and a B8V stellar template. Both fits give similar 
reduced-x^ and are thus equally likely. 
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Table 4. Emission line ratios for the observed and the extinction- 
corrected spectrum of the BCG. The ratio of Ha/H/3 is greater 
than 2.87 due to our assumption that the H/3 emission line is 
depleted due to 2A of stellar absorption. 



4-3.2 Stellar populations 

The BCG has a different spectral shape compared to the 
other cluster galaxies (B2-B6). Whereas the cluster galaxies 
have red B — R colours and high D4000, typical of ellipti- 
cal galaxies, the BCG has blue B — R colours and shallow 
D4000. These spectral features, in addition to the strong 
emission lines, means the BCG hosts a strong source of ul- 
traviolet (UV) emission such as an AGN or a young stellar 
population. 

We determine the stellar composition of the BCG by 
modeling the optical spectral energy distribution (SED). 
Because the BCG is so massive we assume the underlying 
galaxy consists of an old stellar population which has the 
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SED of an elliptical galaxy. The elliptical template from the 
Kinney-Calzetti Atlas (Kinney et al. 1996) is used to model 
this component. 

The blue continuum was modeled by a power-law con- 
tinuum (in which the slope and normalisation were free pa- 
rameters) or early-type stellar templates (05V, 09V, BOV, 
B3V, B5-7V, B8V and AOV) from the Pickles Atlas (Pickles 
1998). The normalisation of the elliptical and stellar tem- 
plates were free parameters. 

The emission lines were masked for the least-squared 
fitting procedure. Since the red part of the spectrum is con- 
taminated by fringing, we only modeled the continuum be- 
tween 3400-5600A and assign constant errors to the entire 
wavelength range of the optical SED. The errors were ad- 
justed to obtain a reduced of approximately x?cd ~1 for 
the best fit model. 

The elliptical plus power-law model resulted in a poor 
fit to the BCG SED (x?ed=13) as the break between 3600- 
4000A was not adequately fit. Thus the blue continuum is 
not due to an AGN component. 

The best fit to the data, with a x?cd— 1-8j resulted from 
the combination of the elliptical template and a B3V stellar 
template. This fit is shown in the top panel of Fig. 7 together 
with the relative contributions of each component to the 
galax;y's luminosity. 

An elliptical plus O-star template is rejected with a 
Xrcd=6.5. The O-star template is similar to the AGN power- 
law SED, and the poor fit results from the same inadequate 
fit to the 3600-4000A break. 

An 05V star emits 5 orders of magnitude more ioniz- 
ing flux than a B3V star, thus to estimate the maximum 
ionizing flux allowed by the optical SED, we flt the optical 
SED with a 3 component model: an elliptical template, plus 
an 05V stellar template, and an additional stellar template 
from the 05V to AOV range. The best flt to the data, with 
a Xred=l-5, results from the elliptical + 05V + B8V model. 
The bottom panel of Fig. 7 shows the model overplotted on 
the optical SED of the BCG, and the relative contributions 
of each component. 

The difference in x?cd between these two models results 
from the poor fit to the Balmer absorption lines. The ob- 
served SED contains strong emission lines so we are unable 
to fit any absorption lines to the data. The B3V stellar tem- 
plate has stronger absorption lines than the 3-component 
model and thus results in a larger x?ed- We therefore cannot 
distinguish which model is more plausible from the Xrod- 

In Table 5 we list the relative flux contributions from 
the 2 and 3 component flts, and the number of ionizing 
stars (steUar types 05V and B3V) that lie in the BCG. The 
monochromatic 4500A luminosities of an 05V and B3V star 
are 1.1 x lO^** and 4.8 x lO^^ergs"^, respectively (Kurucz 
1993). B8V and later stellar types that compose the ellip- 
tical galaxy template emit negligible amounts of ionizing 
photons. 

For both models approximately 60% of the light at 
4500A comes from the young stellar population and only 
40% comes from the elliptical galaxy. However, at longer 
wavelengths the light from the old stellar population dom- 
inates, whilst at 3400A the young population emits ~ 90% 
of the galaxy's luminosity. 

The Ha luminosity resulting from stellar photoioniza- 
tion is calculated from eq. 1 in Allen (1995), using the 



Panagia (1973) ionizing fluxes and assuming unity cover- 
ing fraction. The Ha luminosity expected for the 2- and 

3- component models is 6x10*" and 3.4x10*'^ ergs^^, respec- 
tively. 

The total extinction-corrected Ha luminosity of the 
BCG nebula is 3.4x10*^ ergs"^ therefore between 0.2-100% 
of the Ha emitted by the galaxy is ionized by the young stel- 
lar population. 

The 3-component fit gives an upper limit to the amount 
of ionizing photons emitted by the young stellar population 
because it is forced to contain the largest possible number 
of 05V stars and a covering fraction of unity is asumed. 
Therefore the Ha luminosity from this model is also an up- 
per limit. So whilst it is possible that all the observed Ha 
results from stellar photoionization, there is ample room for 
additional sources to contribute to the ionization of the neb- 
ula. 

4- 3.3 Distribution of the stellar populations and nebula 

Fig. 5a displays an unsharp-mask image of the BCG and B2, 
created by subtracting a smoothed image from the VIMOS 
_R— band image. The galaxy is not smooth, but contains a 
bright filament that traverses the prominent nuclear region 
and extends NE to SW across the galaxy. 

The IFU data is used to visualize the same R-band con- 
tinuum without any emission line contamination from the 
H/3, [Olll] and [Nl] emission lines (Fig. 5b). The centre of the 
BCG is estimated from this image and marked by a cross. 
The red continuum is centrally concentrated in an elliptical 
shape and there are no bright filaments extending southwest 
and northeast. Therefore the large filament must be due to 
the line emission that falls within the R-band passband. The 
nearby galaxy B2 is also visible at position (—20,0). 

Fig. 5c shows the BCG as seen in the shortest wave- 
length emission measured in the IFU spectra (3885-4030 A), 
which in the rest-frame of the cluster is ~3190-3310A, falling 
approximately in the U-band. This wavelength range does 
not include any bright emission lines so the light is emitted 
from the young stellar population. Both NE and SW fila- 
ments are prominent in the blue continuum image, implying 
they are the locations of the recent star formation. 

Fig. 5d maps (in [On]) the emission line nebula that sur- 
rounds the BCG and shows that both the NE and SW fila- 
ments are clearly visible. The brightest region is the galaxy 
nucleus, and the general shape of the nebula follows the blue 
continuum. 

We highlight the differences between the blue and red 
continuum emission with Figs.5e and f which display the 
rest-frame B — R continuum, and the strength of the 4000A 
break, D4000. Fig. 5e shows that the BCG generally has blue 
colours, but the nucleus and SW filament are clearly bluer 
than the rest of the galaxy. We note that the colour varia- 
tions cannot all result from the applied extinction correction. 
The highest Ha/H/3 ratios were observed in the nuclear re- 
gion (see the E{B — V) map in Fig. 1), which translates into 
a large extinction correction. Therefore the enhanced blue 
colours from the nucleus may result from an excessive ex- 
tinction correction, however, the dust is patchy and did not 
extend along the SW filament. 

The variation in D4000 across the BCG is shown in 
Fig. 5f. D4000 is greater than 1.5 in the eastern region of 
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Model 


% E 


% 05V 


% B3V 


% B8V 


no. of ionizing stars 


Ha flux 












xlO'' 


1043ergs~l 


E + B3V 


37 




63 




596 


0.006 


E + 05V H 


h B8V 43 


27 




30 


4.5 


3.4 



Table 5. Stellar populations in the BCG. The optical SED of the BCG is best fit by an old stellar population (modeled by an elliptical 
galaxy template, E), plus a large component of young stars having a B3V SED (E + B3V). In a 3-component fit which is forced to 
include an 05V stellar population, the SED is best fit by a composition of an elliptical galaxy, plus 05V and B8V stars (E + 05V + 
B8V). The percentage of light emitted at 4500A by each component is given in columns % E,% 05V, % B3V and % B8V. The number 
of ionizing stars is the number of 05V or B3V stars, and column 7 lists the Hq flux resulting from photoionization by these stars. The 
E + 05V + B8V model gives an upper limit on the Ho flux produced by stellar photoionization. 



the galaxy, although the region with D4000 > 1.7 at (-20,0) 
is the continuum from the nearby galaxy B2. The central and 
southwest region of the galaxy has a shallow 4000 A break 
with D4000 < 1.3, much lower than observed in elliptical 
galaxies. D4000 is low in the nuclear region (~ 1.1) and 
decreases smoothly down the SW filament of young stars 
to the western tip of the emission line nebula, 42kpc away. 
D4000 is low in these regions, independent of the extinction 
correction, supporting the above finding that the central and 
SW parts of the galaxy host the young stars. 

In summary, the old stellar population of the BCG lies 
within an elliptical shape. Extending beyond this smooth 
profile are two filaments, extending both NE and SW, which 
contain young, blue stars. These filaments are visible in the 
U-hsaid continuum and in line emission. There is a smooth 
decrease in the strength of D4000 for 42 kpc along the SW 
filament implying progressively younger populations along 
the filament, or a larger fraction of younger stars compared 
to the underlying older stellar population. 

4.4 The ionized nebula surrounding the BCG 

4-4- 1 Luminosity 

The flux in the [On] line from the entire ionized nebula is 
9.8 X lO"^** erg s"'^ cm~^ so it is more luminous than any neb- 
ula in the Crawford et al. (1999) sample of BCGs (however 
it should be noted that the Crawford et al. 1999 fluxes lumi- 
nosities are measured within slits and may be incomplete). 
Correcting the spectrum for reddening and extinction boosts 
the [On] luminosity to 4.7 x 10*^ ergs"'^ making it the most 
luminous nebula around a BCG yet observed. 

A Galaxy Evolution Explorer (GALEX) GR4 observa- 
tion gives Mnuv = 18.78 and Mfuv = 18.36 (AB) for 
the BCG. In both NUV and FUV images, the emission is 
spatially extended along the direction of the fllament. The 
observed FUV flux is brighter than expected given the NUV 
flux as dust extinction generally produces NUV>FUV. Since 
the Lya emission line is shifted into the FUV passband, it is 
likely that the observed FUV emission is enhanced because 
of contamination from Lya. 

The UV slope (3, defined as fx = A^, ranges between 
^2.1 and —2.6 for a galaxy which has been continuously 
forming stars for more than a Gyr (Leitherer et al. 1999). 
Almost all other stellar population models result in a higher 
P, therefore we can assume that P > —2.6 in this BCG. 
Thus from the observed NUV flux, we derive an excess flux of 
1.3 X lO"'^^ ergcm"'^ s~^ within the observed FUV passband. 
Correcting for Galactic and intrinsic extinction, we derive a 



lower limit of the Lya emission of 3 x 10~^^ erg cm~'^ s~^ , 
translating into a luminosity of 3.8 x 10'''' ergs"'^. This is 
approximately a factor of 10 greater than the Ha flux, and 
results in a Lya /Ha ratio similar to the theoretical case A 
prediction. 

4-4-2 Mass and volume-filling factor 

The X-ray derived temperature and central electron density 
of the ICM are kT ^ 5.6 keV and rieo = 0.13 cm"^ (see 
Table 2 and Bohringer et al. 2005). In a pure hydrogen, 
completely ionized gas approximation, the ICM pressure is 
therefore ~ 2 x 10~® dyncm"^. Thus, assuming hydrostatic 
equilibrium between the 10, 000 K nebula and the ICM, the 
electron density inside the nebula is no ~ 850 cm~^. 
The mass of the ionized gas is given by 



where rie is the electron density, an^ is the effective re- 
combination coefficient and hvjia is the energy of a photon 
at the frequency of Ha (Osterbrock & Ferland 2006). For 
Uc ~ 850 cm~"^, the equation above yields a mass of approx- 
imately 2.5 x 10'' Mq. 

Assuming this mass is distributed inside a sphere, the 
radius of the sphere is ~ 50 pc. On the other hand, the emis- 
sion in this limited region is spread over a circle of radius 
~ 20 kpc. Assuming again a spherical distribution, the gas 
fills a volume of ~ 3.4 x lO'^'^pc"^. Hence, the filling factor 
is approximately 1.5 x 10~*, a value comparable to those 
found by Hatch et al. (2007) for four BCGs at different red- 
shifts. This low filling fraction suggests either a clumpy or 
filamentary distribution of the optically-emitting gas. 

4-4-3 Source of ionization 

Shocks, conduction from the X-ray emitting ICM and UV 
ionization from hot young stars or the AGN may all play a 
role in heating and ionizing the gas. But the line emission 
may also be heated by secondary electrons from more en- 
ergetic particles or cosmic rays, proposed by Ferland et al. 
(2009) and shown to reproduce the observed line ratios of 
BCG nebulae. Here we study the emission line ratios in order 
to determine the dominant ionization and excitation source 
of the gas. 

Baldwin, Phillips & Terlevich (1981) designed a method 
to discriminate between Hll regions, power-law ionization 
and shock heating. The BPT diagrams (Baldwin et al. 1981) 
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Figure 8. Diagnostic emission-line diagrams showing [Nil] A6584/Ho vs. [Olll] A5007/H/3 and [On] A3727/H/3 vs. [Om] A5007/H/3 in 
logarithmic space. Separation lines used to differentiate between AGN, composite and stellar ionization are from Kewley et al. (2006) 
and Lamarcille et al. (2004), respectively. 



we use are displayed in Fig. 8. The nebula exhibits a vari- 
ety of ratios depending on the location, but no region is 
dominated by AGN heating, and few regions are dominated 
by stellar UV ionization. Rather, these diagnostic diagrams 
show that the gas is ionized by multiple sources. The spa- 
tial distribution of the N[ll]/Ha, 0[ll]/H^ and 0[lll]/H/J are 
shown in Fig. 6. We use these figures to explore whether 
different ionization sources are located in certain regions. 

The [Olll]/H/3 fiux ratio, a good indicator of the ion- 
ization parameter, decreases radially towards the outskirts 
of the nebula, but is slightly offset from the centre of the 
galaxy. The [Olll] may be ionized by the AGN in the nu- 
cleus. This distribution is not shared by any other line ratio, 
but is similar to the distribution of the [On] flux. 

The [Oll]/H/3 ratio is relatively constant around the nu- 
cleus and along the protruding filament. South of this fila- 
ment, the ratio suddenly increases to ^ 0.7, while East from 
the nucleus and towards B2 it decreases below 0.5. These 
low [Oll]/H/3 ratios are consistent with stellar UV photoion- 
ization. If an AGN dominates the ionization in the central 
part of the BCG, not only would we expect higher [Oll]/H/3 
ratios at that position, but also a gradual decrease of this 
line ratio towards the outskirts, where photoionization from 
massive stars should surmount AGN ionization. This is not 
observed, so we rule out this hypothesis. 

The [Oll]/H/3 flux is higher towards the South, suggest- 
ing the existence of a harder, non-stellar excitation source 
in this region. Fig. 9 displays the X-ray surface brightness 
contours from Bohringer et al. (2005) overplotted on top 
of the VIMOS image. Two X-ray maxima are evident; one 
overlapping the BCG center, the other coinciding with the 
region of high [Oll]/H/3. Thus this hard ionization source is 
also visible in X-rays. 

East of the nucleus, the [Nil] /Ha ratio varies rapidly 
between -0.35 and 0.15 (in logarithmic space), which is too 
large to be due to stellar photoionization. Possible heating 
mechanisms include shock heating or power-law ionization. 
To the West, however, this flux ratio is consistently low up to 
a distance of ^ 15 kpc from the centre; then there is a slight 
increase again in the SW tip of the nebula. The region of 




Figure 9. VIMOS R-band image with overplotted contours of 
the X-ray surface brightness from a Chandra image of Bohringer 
et al. (2005). 

low [Nil] /Ha flux ratios coincides with that of a low D4000 
and a blue B-R colour, suggesting the gas in this region 
is ionized by UV radiation from young stars. However the 
BPT diagram (Fig. 8a) indicates that this region is heated 
by more than one source, so an additional source competes 
with the UV from young stars to ionize the gas. 

Fig. 10 shows the three line-ratios as a function of line 
flux. The line ratios saturate at high luminosities, whereas 
fainter regions are characterized by more variability. This 
trend has also been observed in other luminous BCGs by 
Wilman et al. (2006) who postulate that the line ratios sat- 
urate when star formation is proceeding at such a high rate 
that stellar UV dominates the photoionization of the gas. 

In conclusion stellar UV is a plausible ionizing source, 
but it cannot act alone: harder ionizing sources seem to be 
also present. The stellar UV ionization is concentrated along 
the SW filament, whilst the ionizing sources of the nuclear 
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Figure 10. [Oiii] vs. [Oin]/H/3, [On] vs. [Oii]/H/3 and Ha vs. 
[Nii]/Hq. All values on the x-axis are fluxes in units of lO^"*^^ 
erg cm~^ . 



and other regions remain unknown. The ionizing source of 
the nuclear region is of particular importance as this region 
is brightest in both the U-band and line emission. Our data 
point towards ionization by multiple sources, including stel- 
lar UV. 



4-4-4 Gas kinematics 

Fig. 11a presents the line-of-sight kinematics of the nebula, 
as derived from the relativistic Doppler shifts of the [On] 
emission lines. All velocities were computed with respect to 
the median redshift of the [On] A3727 emission {z = 0.2173). 

The galaxy is broken up into distinct regions. The 
galaxy B2 has a gas velocity of — 50kms~^ , and is clearly 
distinct from the red-shifted gas that lies to the East of the 
BOG and in the nuclear region, which has a velocity between 
and 220kms~^ . To the West of the nuclear region the gas 
is blue-shifted up to —220 km s^^ , although at the western- 
most tip of the nebula the gas velocity again changes direc- 
tion and is red-shifted by 200 kms~^ . The central part of the 
SW filament is blue-shifted, whilst the SW tip is redshifted, 
thus the filament must either be stretching or collapsing on 
itself. Similar kinematics can be seen in NGC 1275 in the 
Perseus cluster (Hatch et al. 2006), and in the models of 
Pope et al. (2008). The length of the SW filament is 42 kpc, 
implying a dynamical age of ~ 10* yr, so the line-emitting 
gas must be long-lived. 

The dispersion map of the gas is displayed in Fig. lib. 
All line-widths are cr-widths, and not FWHM (full width at 
half maximum). The cr-widths of the [Olll] and H/3 lines used 
for creating the dispersion map of the gas were corrected for 
instrumental broadening of 375kms~^ as measured from 
the nearby 5755A skyline. White bins in this figure are re- 
gions in which the spectral resolution is too low to determine 
a dispersion. A peak in the dispersion is seen slightly West 
of the nuclear region. The dispersion here is 250-300 km s~^ 
whilst the median velocity dispersion is ~200kms^^ . The 
SW filament has the lowest velocity dispersion. These strong 
internal motions may result from the superposition of dis- 
tinct components along the line of sight. 
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Table 6. Fluxes (J^) are stated in units of 10^^* ergs^^ cm^^, 
FUV flux density is in units of 10~^^ crgs~^ cm~^ Hz~^, lumi- 
nosties (L) are stated in units of 10''^ ergs"'^ and star formation 
rates (SFR) are given in units of yr~^ . 



4.5 Star formation rate of the BCG 

The FUV light from a galaxy is dominated by emission from 
young O and B type stars so we can use the rest-frame FUV 
emission to estimate the star formation rate. This estimate 
of the SFR does not assume that the line emission has been 
ionized by the young stellar population. Since the NUV and 
FUV emission is spatially extended, approximately in the 
same direction as the filament, it is likely to be emitted 
from a spatially extended source, such as stars, rather than 
the central active nucleus. 

Assuming all of the ISOOA flux is emitted from a popu- 
lation of hot young stars, we use the observed NUV emission 
(equivalent to the rest-frame 1870A) to obtain a star forma- 
tion rate. The NUV fiux of 1.1 x lO"^'' ergs"^ cm"^ Hz"^ is 
corrected for Galactic extinction of E{B — V)— 0.108 and 
internal extinction of E{B — V)— 0.211. We use the rela- 
tion between star formation rate and rest-frame FUV flux 
(observed NUV fiux) derived by Salim (2007) assuming a 
Salpeter mass function. 

SFRfvv = 1.08 X 10"^* Lfuv Mq yr~\ (4) 

This equation is calibrated for FUV emission centered 
on 1528A, whilst the rest-frame wavelength of the light ob- 
served through the NUV filter is 1866A. However Salim 
(2007) show that this will have a negligible effect on the con- 
version between FUV luminosity and star formation rate. 
Star formation rates derived from the observed NUV fiux 
are given in Table 6. Kennicutt (1998) derive a relation that 
results in 30% greater star formation rate for a given NUV 
fiux. 

The Hq and [On] luminosity may also be used to es- 
timate the star formation rate in the BCG under the as- 
sumption that the gas is ionized by UV photons from hot 
young stars. Table 6 lists the total fiux, extinction corrected 
fiuxes and luminosities of both [On] and Ha, together with 
the derived star formation rates. Star formation rates (SFR) 
are computed from the [On] and Ha luminosites according 
to the empirical relation of Kewley, Geller & Jansen (2004), 

S-F^wol = 6.58 X 10-*^ M0yr-\ (5) 

and the theoretical equation of Kennicutt (1998), 
SFRl'^gs = 7.8 X 10-"' Lhc Mq yT~\ (6) 
respectively, where ijQjjj and Lua are expressed in units of 
ergs^^. 

The FUV-derived star formation rates (calculated from 
the observed NUV emission) are 2-3 times lower that the 
rates derived from the Ha and [On] emission lines. In sec- 
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Figure 11. Kinematic properties of RXJ1504, from left to right; radial velocity derived from the redshift of the [On] doublet and 
velocity dispersion cr, both in units of kms~^ . The zero velocity point is defined as the redshift of the [Oil] emission line measured in 
the combined spectrum shown in fig 4. 



tion 4.3.2 we show the number of O and B stars present 
in the BCG is barely sufficient to power the Ha emission. 
The discrepancy between the star formation rates derived 
from FUV and emission line luminosity reinforces this con- 
clusion and implies that the [On], Ho? and other emission 
lines are heated by an additional source. Star formation rates 
of BCGs derived from line emission must only be considered 
upper limits of the true star formation rates. 



5 DISCUSSION 

5.1 Brightest cluster galaxy growth 

We present a number of observations which imply the BCG 
is rapidly forming stars. These include extended FUV emis- 
sion, blue galaxy colours atypical of elliptical galaxies, shal- 
low D4000, and some line ratios indicating stellar UV as a 
possible ionization source. The IFU continuum and line ra- 
tio maps indicate this star formation occurs mainly in the 
galaxy core and in a 42 kpc filament stretching from the core 
to the SW. 

The rest-frame FUV emission is likely to be the most 
reliable estimator of the star formation rate, thus setting 
it at ~136 M© yr~^ . This approximately agrees with the 
mass deposition rate derived from the X-ray spectra (~ 
80 M0 yr~^ ). Therefore, the current star formation may be 
fueled by condensing intracluster gas. 

The current star formation rate is large, and together 
with the BCG of Abell 1835, it is the largest star forma- 
tion rate of any low-redshift elliptical galaxy. However, this 
current rate is considerably less than the average rate at 
which stars must have formed over the galaxy's entire life- 
time. The specific star formation rate of the BCG (the ratio 
of current SFR to the mass of galaxy) is ~ 2.3 x 10~^^ yr~^ 
and is therefore similar to many elliptical galaxies which are 
detected in the UV. Had the current rate of star formation 
lasted since z = 1, the galaxy would have increased its mass 
by less than 15%. If this phase of star formation should only 
last a few tens of Myrs, the galaxy would return to its pre- 
starburst colours in a short time after the end of the burst. 
Thus, while the ongoing star formation clearly infiuences the 



galaxy's colours and ionization state of the gas, its impact 
on the mass of the galaxy is not significant. 



5.2 Heating of the nebula 

The optical emission line ratios indicate stellar photoioniza- 
tion dominates in the centre of the galaxy and along the star 
forming filaments. The strong blue continuum and shallow 
D4000 in these regions indicate the presence of ionizing O 
and B stars. 

However, the BPT diagnostic diagrams show stellar 
photoionization cannot be the exclusive ionization source 
and at least one additional mechanism needs to be present. 

Both Ha and [On] line luminosities estimate higher star 
formation rates relative to the FUV derived rate. The num- 
ber of ionizing O and B stars in the BCG, inferred from 
SED fitting, is barely sufficient to produce all the Ha flux 
emitted. This also implies that some of the line emission is 
ionized by an additional source. Thus the Ha or [On] flux 
should only be used to calculate an upper limit of the star 
formation rate. We have also shown that dust extinction 
can be large in BCGs, hence the observed line flux will be 
greatly reduced, and star formation rates derived from single 
emission lines should be considered highly uncertain. 

Line ratio maps are a good means to detect where the 
dominant ionization source changes. The best ratio for this 
work is [On]/H/? because the [Nil] and Hq lines lie in a wave- 
length region that is badly affected by fringing, so the ratio 
of these lines is unreliable. The [On]/H/3 map clearly marks 
the change in dominant ionization source in the South, which 
coincides with a second X-ray peak. None of the line ratio 
maps show a correlation with the map of recent star forma- 
tion (Fig. 5c), providing further evidence that stellar pho- 
toionization is not the only ionizing source at work. This 
result does not come as a surprise and it has also been ob- 
served in other BCGs (e.g. Sabra et al. 2000; Hatch et al. 
2007). Ferland et al. (2009) have proposed an alternative 
particle heating method which reproduces the line ratios. 
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5.3 AGN Feedback 

The mass deposition rate allowed by the X-ray spectra is 
less than 10 percent of the mass deposition rate of 1400-1900 
Mq yr~^ implied by the total X-ray luminosity emitted from 
the cluster core. Since we have an indication that 10 percent 
of the ICM is cooling and forming stars, it follows that a 
heating source must supply 90 percent of the luminosity 
emitted from within the cooling radius, which is 2.7 x 10^''' 
erg s~^ . If this is supplied by AGN feedback from the central 
black hole, the mass accretion rate must be 0.5 Mq yr~^ , 
assuming at least a 10 percent output efficiency. 

The K—hand magnitude of a galaxy correlates with the 
mass of the supermassive black hole (SMBH) at its centre 
(e.g. Marconi & Hunt 2003; Graham 2007): 

log (Mbh/ Mq) = -0.33 [Mk + 24] + 8.33, (7) 

where Mbh is the mass of the black hole and Mk is the K- 
band magnitude of the BCG. With Mk = -27.0 (2MASS), 
it follows the BCG has a SMBH of mass ~ 2 x 10^ Mq. 

To estimate the power supplied by an AGN, we first 
consider Bondi accretion from the hot gas (Bondi 1952). The 
Bondi accretion radius is 12 pc. As for all SMBH except 
perhaps the one in M87, this is far beyond the resolution 
limit of our instruments and we can only speculate about 
the physical conditions close to the black hole. The black 
hole will accrete mass at a rate Meondi ~ TrAcspr^, where 
A depends on the adiabatic index (here "/ — 5/3, A = 0.25; 
Bondi 1952), Cs is the sound speed, p is the gas density 
and va is the accretion radius. For the central parameters 
determined from the X-ray data, namely kT = 5.6 keV and 
rie = 0.13 cm"^, we obtain Meondi « 5 x 10"" M© yr"^ . 
However, this temperature and density represent an average 
over a large extraction region used for the X-ray spectral 
analysis, therefore the true temperature and density near 
the black hole could be very different. If we assume that 
the thermal gas pressure at the centre is no less than the 
thermal pressure of the inner region measured spectrally to 
be at a temperature of 5.6 keV, then p scales as 1/T (in 
practice the extra weight will make it go steeper than this). 
Therefore, A^Bondi oc 1/r^'^. So to increase MBondi by 10^, 
to 0.5 Moyr"^ , requires only a drop in temperature by a 
factor of 16 which means to 350 eV. Many nearby cool core 
clusters show a cool component at ~ 0.5 keV (e.g., Sanders 
et al. 2010) and we cannot rule out such a low temperature 
in the centre of RXCJ1504. 

One could also try to infer the AGN power output over 
the past 10** years from the radio power at 1.4 GHz using 
the relation of Birzan et al. (2008): 

log Pjet = (0.35) log Pi4(,0 + 1.9. (8) 

In the relation above, P1400 represents the bolometric radio 
luminosity for the total source at 1.4 GHz in units of 10^* 
WHz^^, and Pjet is the jet/cavity power in units of lO'*^ 
ergs~^. 

The BCG contains a radio source with a brightness of 
62 mjy at 1.4 GHz (Bauer et al. 2000). At a luminosity dis- 
tance of 1075.5 Mpc, this implies a radio power of 8.6 x 10^* 
WHz~^, thus a jet power of 1.5 x 10"''* ergs~^. The corre- 
sponding black hole accretion rate at 10 percent efficiency is 
0.03 M0yr~^ and an order of magnitude smaller than the 
power required to heat the ICM. However, the relation found 



by Birzan has a scatter of 0.85 dex (factor of ~ 7) at 1400 
MHz, so that our object lies a little over la away from the 
mean (~ lAa in log). 

In conclusion, we cannot rule out that heating is sup- 
plied by a central AGN that accretes with a rate of ~ 0.5 
M© yr-i . 

If the accretion rate has been as high as 
0.5 M0 yr~^ since redshift 1, then the black hole has 
grown by an amount of ~ 3 x 10^ Mq, which is larger 
than the total estimated mass of ~ 2 X lO'' Mr;i . This 
implies that the current accretion rate is atypically high, 
which suggests that the present phase in the feedback 
cycle with the high star-formation and accretion rate must 
be relatively short lived. Furthermore, BCGs may harbor 
supermassive black holes with masses a few times higher 
than those predicted by the Mbh — Mk relation of Graham 
(2007), as is the case in Abell 1836 (Dalla Bonta et al. 
2009). 

5.4 Turbulent energy of the ICM 

Assuming the velocity dispersions and the line-of-sight ve- 
locities of the cold gas are solely due to turbulence gives an 
upper limit on the turbulent velocities in the nebula. Despite 
the large density contrast between the ICM and the optical 
line-emitting gas, the low filling factor suggests that the neb- 
ula gas may be blown about by the hot gas and thus serves 
as a good tracer of ICM motions (Fabian et al. 2003). The 
speed of sound in the hot central plasma at a temperature of 
5.6 keV is roughly 1190 kms~^ and the characteristic veloc- 
ity of isotropic turbulence is vturh = •\/2 (u^ -I- Va), where 
?;a is the dispersion of the velocity shear of the cold gas, 
measured from the a-width of the line-of-sight velocity his- 
togram, wa ~ 117 kms^^ , and Va is the average velocity 
dispersion, ~ 200 kms~^ . iiturb provides an upper limit 
on the velocity dispersion of the hot gas. vturh ~ 325 km s~^ 
and therefore the turbulence in the central ICM has a typ- 
ical Mach number of less than 0.27. We thus obtain an up- 
per limit on the ratio of turbulent to thermal energy in the 
ICM, given by (7/2) M^ with 7 = 5/3 the adiabatic index 
of monatomic ideal gas, of roughly 6%, in agreement with 
upper limits derived from other methods (Churazov et al. 
2008; Werner et al. 2009). 



6 CONCLUSIONS 

RXCJ1504. 1-0248 is one of the most extreme cooling flows 
ever discovered. The observations presented here offer the 
first multi-wavelength view of the cluster core, examining 
both the hot and warm gas as well as the stars of the BCG. 

The core of RXCJ1504.1-0248 has an X-ray luminosity 
of at least 3 x 10^^ erg s~ ^ , approximately 10 percent of which 
is supplied by gas cooling and forming stars. Accretion onto 
a supermassive black hole of 2 x 10® Mq is the likely mecha- 
nism needed to account for most of the remaining 2.7 x 10*^ 
ergs~^ heating. The corresponding Bondi accretion rate at 
10 percent radiative efficiency is 0.5 Mq yr~^ . Using the 
correlation between 1.4 GHz radio luminosity and jet power 
(Birzan et al. 2008), we infer a jet power of 1.5 x 10** ergs~^, 
still insufficient to heat the ICM. 

The intracluster medium mass deposition rate is in 
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agreement with the star formation rate of the BCG derived 
through FUV observations (~140 Mq yr"^ ). The IFU im- 
ages show that most of the star formation occurs in the 
galaxy core and in a 42 kpc-long filament that stretches out 
from the core to the SW. These regions have blue continuum 
colours, and shallow 4000 A breaks consistent with recent 
star formation. Additionally, the [On]/H/3 line ratio indi- 
cates stellar photoionization is the dominant source of ion- 
ization along the filament. The X-ray emission also extends 
along the filament indicating a direct connection between 
the intracluster medium and the recent star formation. 

The line-emitting nebula that surrounds the BCG is the 
most luminous observed to date, with extinction-corrected 
[On] and Ha luminosities of 4.7 x IC*^ ergs~^ and 3.4 x 
10''"^ ergs~^ respectively. The number of ionizing O and B 
type stars in the BCG is barely sufficient to power this bright 
nebula, and the BPT diagrams and line ratio maps show 
that an additional heating mechanism is required. Because 
this heating source produces Ha emission, BCG star forma- 
tion rates that are derived from Ha luminosities should be 
considered upper limits. 

The nebula kinematics are ordered and the nebula com- 
prises of a number of kinematically distinct regions. The ve- 
locities are generally low (< 250 kms~^), and there is no 
evidence for rotation or free-falling gas infall. The 42 kpc 
filament has a velocity shear of ~400kms~^ resulting in 
a dynamical time of 10* years, thus the nebula must be 
long-lived. The average velocity dispersion of the nebula is 
Vcr ~ 200 kms~^ and the dispersion of the velocity shear of 
the cold gas is «a ~ 117 kms^^ . These set an upper limit of 
~ 325 kms~^ on the mean turbulent velocity in the ICM, 
meaning that the ratio of turbulent to thermal energy of the 
intracluster medium is limited to less than 6 percent. 

Our results suggest mass transfers from the intracluster 
medium to the brightest cluster galaxy at a rate of approxi- 
mately lO'^ Mq yr~^ . The cooling of the intracluster medium 
results in a gas reservoir in the brightest cluster galaxy. The 
rate at which gas condenses into this reservoir is similar to 
the rate at which this reservoir forms stars, so the reservoir 
should be approximately stationary. Our observations can- 
not reveal the coolest gas components of this reservoir, but 
we find that the lO" K gas is partly heated by a source other 
than the forming stars. 
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